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ABSTRACT 


The  developing  flow  in  the  hydrodynamic  entrance  region  of  an  equi- 
lateral triangular  channel  was  investigated,,   Using  a  fully  developed 
velocity  profile  from  Knudsen  and  Katz  (3)  as  a  boundary  condition8  the 
equations  were  solved  numerically  employing  the  method  of  Chorin  (l)e 

The  resulting  velocity  profiles  and  pressure  drop  were  employed 
in  calculating  the  local  friction  factor  in  the  hydrodynamic  entrance 
region . 
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I o   INTRODUCTION 

The  flow  characteristics  of  non-circular  channels  have  become 
of  increasing  importance  in  the  analysis  of  compact  heat  exchangers 0 
Nevertheless,  heat  transfer  analyses  of  the  triangular  entrance  region 
have  been  limited  to  approximations  which  postulate  fully  established 
velocity  profiles  (Graetz  Approximation) »   Such  solutions  have  been 
obtained  by  Lu  and  Miller  (4)  and  Tao  (8)  for  equilateral  triangles 0 

Using  Sparrow's  (7)  results  for  the  fully  established  velocity 
profiles  in  an  isosceles  triangular  duct,  McComas  (5)  has  estimated 
values  of  entrance  lengths  and  friction  factors  for  such  flows o 

In  order  to  complete  the  analysis  of  heat  transfer  in  the  entrance 
region  of  triangular  ducts,  it  is  necessary  to  obtain  a  hydrodynamic 
solution  for  the  developing  velocity  profiles  in  that  region.   It  is 
therefore  the  purpose  of  the  present  work  to  provide  a  solution  of 
the  equations  of  motion  in  the  hydrodynamic  entrance  region  of  an 
equilateral  triangular  duct* 


11 


II 9   ANALYSIS 

In  an  attempt  to  find  a  convergent  solution  to  the  equations  of 
motion,  several  explicit  methods  were  tried  without  success s   Subse- 
quently a  relaxation  method  employing  an  artificial  density  change^ 
suggested  by  Chorin  (1),  was  adoptedo   This  method  essentially  alters 
the  solution  of  the  equations  of  motion  from  an  initial  value  problem 
to  a  boundary  value  problem,  necessitating  the  use  of  entrance  length 
and  the  fully  developed  velocity  profile  as  boundary  conditions , 
Ae   Gove  rn  in g  E  q  ua  t  i ons 

The  governing  equations  for  the  flow  of  an  incompressible  viscous 
fluid  (Navier  Stokes  Equations)  may  be  written: 

x  -  Momentum: 

^      c>X       ^  <*<!       ?  ciX  (1) 

y  -  Momentum: 

«ft    &%  ^      ^  ?  ^  (2) 

z  -  Momentum: 

Jt     <*%     ^     «a^     ?  <^  (3) 


Continuity ! 

^       *i         ^  (4) 


T    —    -+-    —  tj 
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Introducing  the  dimensionless  variables, 


v-Z 

V-  £        w-  £ 

X=J 

Y«J        ?-£ 

P  -  -p  -4- 

'■'J 

(5) 


equations   1-4  become: 
it'       "  Jl       •     -M  J-  J  &  (6) 


if 


(7) 


(8) 


&        ^Y  3?  (9) 


B.   Boundary  Conditions 

A  condition  of  uniform  velocity  is  assumed  at  the  entrance  and 
no-slip  conditions  are  imposed  at  the  walls;  thus; 

1«   at  the  entrances* 

B(o.Y.?.0-  1 

V (OJJ.O-  0 

WIO.Y,?,*)  ■  0 

2»   at  the  walls s 


*The  prime  has  been  dropped  and  dimensionless  variables  are 
understoode 
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C„   Artificial  Density 

Because  the  problem  is  one  of  incompressible  flow,  no  density 
change  can  take  place 0  One  can  introduce  an  auxiliary  continuity 
equation,      3f  .  ^U   ,  W  t  IV 

5t*3x     «     *  (9a) 

in  which  changes  of  the  pseudo  density  given  by 

P=f/S  do) 

can  be  used  to  test  for  convergence  at  each  iteration  process 6 
Thus,  we  have  introduced  a  set  of  auxiliary  parameters,  ^  and  £    , 
which  are  analogous  to  an  artificial  density  and  an  artificial  com- 
pressibility,, 

By  introducing  the  artificial  equation  of  state  (equation  10) 
and  using  t'  as  an  auxiliary  parameter  analogous  to  time  in  compressible 
flow,  the  initial  value  problem  can  be  altered  to  a  boundary  value 
problem0   Such  a  change,  however,  requires  the  use  of  a  fully  estab- 
lished solution  as  an  additional  boundary  condition,   An  exact  solution 
for  the  steady,  fully  established  profile  has  been  obtained  by  Knudsen 
and  Katz  (3) s 

D-Cr.i)--2fR[i(7»*f)-i(?V3",j:-^j 

wherei  7  =  Y"jyT  ands  f~l.  "  33  0   The  friction  factor  f  for  fully 
developed  flow  is  given  by  Kays  (2)  as  130330   This  number  serves  to 
determine  the  value  of  the  fully  established  pressure  gradient „ 
Do   Finite  Difference  Equations 

Equations  6-9  were  put  into  finite  difference  form  using  central 
differences,  except  for  the  second  derivatives  which  were  differenced 
using  the  Dufort-Frankel  pattern0   This  latter  expression  is  of  the  form; 
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n  kl  Nl  |sl-t-l        N-' 

^  ^  (12) 

where  UN(X+)  denotes  U(N,X+AX,Y,Z),  UN+1  is  U(N+l,XtY,Z)  and  N 
denotes  the  Nth  relaxation •   The  resulting  equations  were  solved  for 
the  N+l  terms  and  the  resulting  computing  equations  are  of  the  form: 


-  Rf£  [uS.>t&>  - tJ5-)U<?-.J  -  Rg  [uwV^i  -  Uw->V«->] 

-  Rf*  [u^)W*->  -  1&-.V&  ■>]  +  %£  [oV>  *  U «->  -  tf  "J 


-  Ra  jV«.>w?V\r«>We->]  +  ^  [v<k+V«->- V"'  J 


(14) 


wm,-[i/(i^^^)]*[-KUf^>-^] 

+  fc£  [w-«t)  *¥&■)  -VJ  *■  ^  [WcS^w&-,-  w-'j]  ( 


15) 


J 


t-Ui 


(16) 
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Eo   Grid  Size  and  Stability 

The  geometry  of  an  equilateral  triangle  presents  a  problem  in 
selecting  appropriate  values  of  AY  and  AZ„   The  requirement  that  the 
network  intersections  fall  on  a  boundary  was  considered  essential  for 
numerical  calculation^  dictating  the  relationship  between  aY  and  AZ 
such  that  aY  =V3*  &i/Z  . 

The  value  of  AX,  while  arbitrary,  is  constrained  by  the  fact 
that  the  number  of  X-wise  stations  must  be  chosen  such  that  the  number 
of  stations  times  the  value  of  aX  be  equal  to  the  predicted  entrance 
length,  or  AX  x  I=Xe,  where  I  is  some  integer,.   Finally,  because  of 
the  limited  storage  available  in  a  computer,  the  number  of  grid  points 
is  limicedo   For  the  purposes  of  this  investigation,  the  values  adopted 
were  aX=«C0159,  aY=0o 17321,  and   AZ-0o3o   This  grid  resulted  in 
25  stations  in  the  X-direction,  11  in  the  Y  and  6  in  the  Z,  as  shown 
in  Figure  10 

Having  selected  the  grid  size,  the  stability  of  equations  13-16 
may  be  determined  using  the  results  of  Chorin  (1),  who  points  out  that 
the  introduction  of  an  artificial  density  results  in  an  artificial 
speed  of  sound  of  the  forms 

C%*    '/§  (17, 

A  Mach  number  relative  to  the  speed  of  sound  is  given  by, 

M=  5(tC*Vl*W^f  (18) 

For  convergence  it  is  necessary  that  the  Mach  number  (M)  be  less  than 
one6  With  this  condition  satisfied,  the  system  of  equations  is  stable 
when  2 


.       2.  /    Miwio*UM      \V/2. 


(19) 
16 


In  the  present  work  these  values  were  taken  to  be  M=00975,   5=9,38025, 
At=0, 00035,  and  R=lo0„ 
F.   Friction  Factor 

The  wall  shear  stress  may  be  related  to  the  local  pressure  gradient, 
which  is  of  importance  to  engineers  in  predicting  pressure  drope   The 
average  stress  at  any  section  is  given  by  the  Newtonian  shearing  law, 


T.vft). 


(20) 
where  (j^j    represents  the  average  velocity  gradient  evaluated  at 


'8 

the  wall. 


Using  equation  20,  a  dimensionless  friction  factor  can  be  defined, 

2?'*  (21) 


9=i 


in  which  f  is  the  usual  flow  friction  factor,, 
Ge   Numerical  Solution 

m  m  win    HiiWi  i  — MMMMM  m  Ml  r    ■ 

In  the  present  analysis,  using  the  values  of  &X,   aY,  aZ,   ^t 
and  §  given  above,  equations  13-16  were  found  to  converge  to  a  steady 
solution  after  2,000  iterations,,  While  Chorin  points  out  that  initial 
values  of  variables  may  be  set  at  zero  throughout 0  the  convergence  can 
be  hastened  by  using  a  good  "guess"  at  the  values  to  be  calculated,, 
In  the  present  case  the  fully  established  velocity  values  were  used 
throughout  as  an  initial  assumption 

A  copy  of  the  program  written  in  Fortran  IV  language  is  given  in 
the  AppendiXo 
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IIIo   RESULTS 

The  dimension less  centerplane  velocity  profiles  are  given  in 
Figure  2  for  five  values  of  X/Xe  (062,  0oA,  006e  008,  and  lc,0)c   Figures 
3  through  7  show  the  centerplane  velocity  profiles  and  velocity  contours 
for  X/Xe=0629  06At  006D  0680  and  lo0  respectively 0   The  secondary  ve- 
locities, V  and  W,  are  shown  in  Figure  8  at  a  station  about  midway  down 
the  entrance  region  at  X/Xe=0e60   The  value  of  friction  factor  at  each 
x-wise  station  is  plotted  in  Figure  9  in  the  form  of  f/fe  versus  X/Xe> 
where  the  e  subscript  denotes  the  fully  established  valuesc 

The  dashed  line  on  the  friction  factor  curve  was  computed  from 
the  local  skin  friction  predicted  from  a  boundary  layer  analysiso   The 
value  of  skin  friction  must  of  course  be  unbounded  at  X=0  because  of 
the  discontinuity  in  velocity  at  the  lip  of  the  duct0   For  large  values 
of  X/Xg,  f/fe  must  approach  unity  asymptotically  as  the  results  indi- 
cate* 
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IV,   CONCLUSIONS 

From  the  results  of  this  study,  it  may  be  concluded  that; 

1,  The  relaxation  procedure  used  in  this  analysis  was  found  to 
be  convergent  subject  to  certain  stability  parameters  as  Chorin  (1) 
points  out <>   The  results,  however,  remain  unchecked  experimentally • 

2.  The  use  of  the  fully  established  values  as  exit  boundary 
conditions  forces  the  relaxation  values  to  converge  on  these  values 
in  the  limit  as  the  number  of  iterations  increase,,   Therefore,  the 
values  of  the  computed  velocity  profiles  can  be  only  as  accurate  as 
the  exit  boundary  conditions ,   Due  to  the  lack  of  analytic  and  experi- 
mental results  in  the  area,  there  was  little  choice  in  using  the 
results  of  Knudsen  and  Katz  (3)  for  the  fully  established  velocity 
profile  and  Kays  (5)  for  the  fully  established  pressure  gradient, 

3.  In  general,  the  grid  sizes  used  should  be  reduced  until  the 
number  of  points  reaches  the  limit  of  computer  storage.  This  serves 
to  maximize  the  accuracy  of  the  results  consistent  with  the  accuracy 
of  the  fully  established  profiles  used0 

4,  It  seemed  intuitively  obvious  when  first  using  Chorin' s  method 
that  because  it  represents  a  convergent  numerical  solution  to  the 
finite  difference  form  of  the  Navier  Stokes  and  continuity  equations, 
using  an  entrance  length  greater  than  a  predicted  value  would  lead  to 
the  fully  established  profile  upstream  of  the  boundary  condition. 
This,  however,  was  not  the  case.   Thus,  the  necessity  of  using  the 
predicted  value  of  Xe  reduces  the  generality  of  this  method  in  the 
solution  of  physical  problems.   The  use  of  a  Reynolds  number  of  unity 
may  have  been  the  cause  of  this,  however  it  would  seem  evident  that  a 
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complete  analysis  of  the  equations  must  be  made  numerically  to  find 
the  reason  for  this  limitation 

50   The  Reynolds  number  dependency  of  the  governing  equations  can 
be  extracted  by  using  non-dimensional  space  parameters  of  the  form 
X=xR/d0   An  analysis  of  the  resulting  equations  would  seem  beneficial 
as  the  results  would  be  completely  general  and  would  represent  "uni- 
versal" profiles o 
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FIGURE  I. 
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AND  MESH  POINTS. 
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FIGURE   2. 

DIMENSIONLESS  VELOCITY   PROFILES 
FOR  VALUES   OF    X/Xt  . 
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FIGURE  3. 

VELOCITY  PROFILE  AND  CONTOURS 
FOR    X/Xf=0.2 
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FIGURE  4. 

VELOCITY  PROFILE  AND  CONTOURS 
FOR    X/Xf*0.4 
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FIGURE  5. 

VELOCITY   PROFILE  AND  CONTOURS 
FOR   X/X#  =  0.5 
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FIGURE  6. 


VELOCITY   PROFILE  AND  CONTOURS 
FOR    X/X  =0.8 
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FIGURE  7. 

VELOCITY  PROFILE  AND  CONTOURS 
FOR    X/X    =1.0 
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FIGURE  8. 


SECOND  VELOCITY  PROFILES 
FOR    X/Xe*0.6 
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Friction  Foctor  Calculoted  from  boundory  layer  theory  for 
small   X/X. 
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Result*  of  the  Numerical  Analysis 
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FIGURE  9. 
FRICTION  FACTOR   RATIO  vs.    X/Xe 
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APPENDIX 
SAMPLE  FORTRAN  IV  PROGRAM 
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